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Abstract 

The Dirac monopole string is specified for anti de Sitter cosmological model. Dirac equa- 
tion for spin 1/2 particle in presence of this monopole has been examined on the background 
of anti de Sitter space-time in static coordinates. Instead of spinor monopole harmonics, the 
^vq i technique of Wigner /^-functions is used. After separation of the variables radial equations 

have been solved exactly in terms of hypergeometric functions. The complete set of spinor 
wave solutions ^ e ,j, m ,x(t, r, 9, <fi) has been constructed, the most attention is given to treating 
\^ | the states of minimal values for total moment quantum number j m in- At all values of j, the 

i energy spectrum is discrete. 

PACS numbers: ll.10.Cd, 04.20.Gz 

1 Introduction 

De Sitter and anti de Sitter geometrical models are given steady attention in the context of 
developing quantum theory in a curved space-time - for instance, see in [1]. In particular, the 
■ problem of description of the particles with different spins on these curved backgrounds has a 

long history - see [2-34]. Here we will be interested mostly in treating the Dirac equation in 
de Sitter model. In the present paper, the influence of the Dirac monopole string on the spin 
1/2 particle in the anti de Sitter cosmological model is investigated. Instead of spinor monopole 
harmonics, the technique of Wigner L>-functions is used. After separation of the variables radial 
equation have been solved exactly in terms of hypergeometric functions. The complete set 
of spinor wave solutions ^e,j,m,\(t, r ,@i4>) nas been constructed. Special attention is given to 
treating the states of minimal values for total moment quantum number j m i n , these states turn 
to be much more complicated than in the flat Minkowski space. At all values of j, the energy 
^ ' spectrum is discrete. 

W. 

2 Dirac particle in the anti de Sitter space 

The Dirac equation (the notation according to [39] is used) 
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in static coordinates and tetrad of the anti de Sitter space-time 
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From dH), with the substitution ^(j;) = r -1 ^ 1 / 4 ip(x), we get 
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Below the spinor basis will be used 

I 

1 



-CTj 
(Tj 



12 



(T3 

a 3 



3 Separation of the variables 
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Let us start with the monopole Abelian potential in the Schwinger's form |35J in flat Minkowski 
space 



A a (x) = (A», A*)=l0,g 



[r x n) (r n) 



r (r 2 — (r n) 2 ) / 

Specifying n = (0, 0, 1) and translating the A a to the spherical coordinates, we get 
A = , A r = , Ag = , As = g cos . 
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It is easily verified that this potential^ obeys Maxwell equations in anti de Sitter space 
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Correspondingly, the Dirac equation in this electromagnetic field takes the form 
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i/>(x) = , 



where 



2 id^ + (ia — k) cos 9 
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and k = eg /he. As readily verified, the wave operator in Q commutes with the following three 
ones 
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J 3 - l 3 , 
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which obey the su(2) Lie algebra. Clearly, this monopole situation come entirely under the Schrodinger 
[36j and Pauli |37j approach (detailed treatment of the method was given in |40j). Correspond- 
ingly to diagonalizing the J% and J 3 , the function ip is to be taken as (D a = D J _ m a ((f), 9, 0) 
stands for Wigner functions 
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Further, with the he;p of recursive relations |38| 

deD k+1 / 2 = aD k _ 1 / 2 - bD k+3 / 2 , doD k _ 1 / 2 = cD k _ 3 j 2 ~ aD k+1 / 2 
-m-(k + 1/2) cos 9 } D k+1/2 = (-aD k _ 1/2 - bD k+3/2 ) 
-m-(k- 1/2) cos 9 } D k _ 1/2 = (-cD k _ 3/2 - aD k+1/2 ) 



sin 
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a = -V(j + l/2) 2 -k 2 
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we find how the ^ acts on ij) 
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hereafter the factor w (j + 1/2) 2 — k 2 will be denoted by v. For the fi(r), the radial system 
derived is 



r 
v 



dr 



■4/1 + » v^- A + i- h 

-7=h- % V$— / 2 - * - /l 

V<3? c« r 



M /i = , 
M / 2 = , 
Mf 3 = 0, 
M U = . 
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Else one operator can be diagonalized together with idt,J^,J^: namely, a generalized Dirac 
operator 



£fc = _, 7 7 3 S ^. (15) 

From the equation K k i[) e j m = A V'ejm we find two possible eigenvalues and restrictions on fi(r) 

(16) 



h = *h, h = Sf 2 , \ = —6 + 1/2) 2 - k 2 
Correspondingly, the system (fT4"|) reduces to 



where 



/i + / 2 A - / 2 



(17) 



/ 



5 = 



V2 ' * V2i 
It is known that quantization of k = eg /he and j is given by 

eg/he = ±1/2, ±1, ±3/2,...; 
j=\k\ -1/2, | fc | +1/2,| k | +3/2, . . . (18) 

The case of minimal value jmin =| k | —1/2 must be treated separately in a special way. For 
example, let k = +1/2, then to the minimal value j = there corresponds the wave function in 
terms of only (t, r)-dependent quantities 



^=+1/2^) = e ~ 
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At k = —1/2, we have 
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Thus, if k = ±1/2, then to the minimal values j m i n there correspond the function substitutions 
which do not depend at all on the angular variables (9, <p); at this point there exists some formal 
analogy between these electron-monopole states and S-states (with I = 0) for a boson field of 
spin zero: $z=o = 3>(r, t). However, it would be unwise to attach too much significance to this 
formal similarity because that (#, (/^-independence of (e — g)-states is not a fact invariant under 
tetrad gauge transformations. 

In contrast, the relation below (let k = +1/2) 



,+1/2 ,(j=0) 
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n =+1/2 ( x ) =7 cot 9 (iu - r k=+1/2 

is invariant under arbitrary tetrad gauge transformations. Correspondingly, the matter equation 
([9]) takes on the form 



a/c£> at or J 



(22) 



It is readily verified that both (fT9|) and ([20]) representations are directly extended to (e — g) 
states with j = j m ; n at all the other k = ±1, ±3/2, . . . Indeed, 

/i(r) D k _ 1/2 
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and the relation Hg^ipj^ = still holds. For instance, let us consider in more detail the case 
of positive k. Using the recursive relations 

deDk-1/2 = \^2k - 1 D k _ 3/2 , 



sin" 1 9 [ -m - (k - 1/2) cos 9 } D k _ x 
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in a sequence, the identity Eg^ ipj min = holds. The case of negative k can be considered in the 
same way. Thus, at every k, the j m i n -state equation has the same unique form 

('^l + ^s- M )*"' = 0; (25) 

which leads to the same unique radial system 
k = +1/2, +1,... 

-^=h-i h ~ M ft = , 

-J= / 1 + ^i-/ 1 -M/ 3 = 0; (26) 

fc = -1/2,-1,... 

/ 4 + % V*-^ h ~ M h = , 

^= / 2 - * v 7 *^ / a - M / 4 = . (27) 

In the limit of flat space-time, these equations are equivalent respectively to 
k = +1/2, +1,... 

(^ + ^ 2 -m 2 )h = 0, fs = ^(e + 4)h; (28) 



\dr 2 J ' m \ dr 

■1/2,-1,... 

^ + 6 2 -m 2 V4 = 0, / 2 = lf e + 4) /4- (29) 



dr 2 J ' m \ dr 

These equation both lead us to the functions / = exp(±\/m 2 — e 2 r). In particular, at e < m, 
we have a solution 

exp ( — \J m 2 — e 2 r ) , (30) 
which seems to be appropriate to describe bound states in the electron-monopole system. 

4 Solution of the radial equations 

Let us turn back to the system (|17[) and (for definiteness) consider equations at 5 = +1 (formally 
the second case 5 = — 1 corresponds to the change M =^> —M) 

(v^-f + - ) / + (4= + M) g = , 
dr r -v/cj> 



(v^-f- - - ) g - (-^ - M) / = . (31) 
dr r J(& 
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Here we see additional singularities at the points 

e + V¥ M = or e - V$ M = 
For instance, the equation for f(r) has the form 
2r Mr 
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/ = o. 



However, there exists possibility to move these singularities away through a special transfor- 
mation of the functions f(r),g(r) |24j . To this end, let us introduce a new variable r = sinh p, 
eqs. (f3"T|) look simpler 



d ^ 
dp sinh p 

, d 



)/ + 
)*-( 



cosh p 
e 



dp sinh p cosh p 
Summing and subtracting two last equations, we get 

€ 



+ M)g = 0, 
- M)f = . 



Uf + g) + -^-(f-g)- , 

a/9 sinn p cosn p 



dp ^ ^ sinh p 



(/ + <?) + 



cosh p 



(/ - 5 ) + M(f + g)=0, 
(f + g )-M(f-g) = 



Introducing two new functions 

f + 9 = e-t>l\F + C) - 

or in matrix form 



= e 


p/2( 


G 




H 





/-<? 



hp/2 



(F - G) 



cosh p/2 — sinh p /2 
— sinh p/2 cosh p/2 



where (see definition of the variable z below) 



cosh — 
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one transforms (l33l) into 
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sinh — 
2 



(32) 



(33) 



(34) 



(35) 



(36) 



— e' p/2 
dp 



(F + G) + 



sinh p 



(F-G) 



coshp 



_ G ) + Me- p/2 (F + G) = Q, 



^ e +p/2 (F _ G) + _ 
ap sinn p 



e -p/2( F + G ) 



+-^— e-^ 2 (F + G) - Me+^ 2 (F - G) = 
cos p 
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or 



|(F + G )-I(F + G ) + -^- 
e 



(cosh p + sinh p) (F — G) 



cosh p 



(cosh p + sinh p) (F - G) + M(F + G) = , 



-y-(-F - G) + ^(F - G) + — ^-(coshp - sinhp)(F + G) 
ap 2 sinh p 



+ - 



-(coshp - sinhp)(F + G) - M(F - G) = . 



cosh p 

Now summing and subtracting two last equations, we obtain 

d ^ cosh p sinh p ^ ^ ^ ^ ^ ) G 

<ip sinh p cosh p 2 

. d cosh p sinh p . _ , _ _ 1 . _ 

— - * -r-r^ + e — f-) G + -e + M + - - F = 

dp sinh p cosh p 2 
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Let us translate eqs. (|3Tj) to the variable 2: 



r 2 = sinh 2 p 



d / ; 



, ~, r d \f\ — z yj—z 

2y/-z(l-z) — + u V - - e ; 
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1 



{ 2 ^^Tz 



+(e + M-u--)G = 0, 



v , — +e 



G 



+(-e + M + v- -)F = . 

From (138p it follow two 2-nd order differential equations for F and G respectively 
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dG 
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G = 0. 



4(1 - z) 4z 
With the use of substitutions 

F = z A {\ - z) B F{z) , G = z K {l-z) L G(z), 



eqs. take the form 



z(l-z) 
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d 2 F 
dz 2 



2A + -- (2A + 2B + l)z 



dF 

dz 



(A + Bf 



e(e- 1) + 2B{1 - 2B) 



v(y + l) - 2A(2A - 1) 

Iz 



z(l - z) + 



dz 2 



2K + -- (2K + 2L + l)z 



dG 

dz 
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(K + L) 



e(e + 1) + 2L(1 - 2L) 



4(1 - z) 

v{v - 1) - 2K (2K - 1) 
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First let us consider eq. PU|) ; at A and B taken accordingly 

l + v 



G = 0. 



it becomes simpler 



.4 



v 
2 ' 



e 1 - e 
2' 2 



2A + - - (2A + 2B + l)z 



d£ 

dz 



+ 



- I M - - 



{A + Bf 



f = 0, 



which is of hypergeometric type with parameters 

_ M 1 
a ~ T ~ 4 

To construct functions appropriate to describe bound states we must choose 



+ A + B 



Ml, A 
— + - + A + B, c = 2A + 12. 
2 4 



A 



l + v 



>0, 



B 



1 



<0, c = v + 3/2; 



2 2 
polynomial solutions will arise with the quantization rule imposed 



a = —n , 
b = - n - M -1/2 



e n = M + 2n + v+- 



c = v + 3/2 . 

Now let us turn to eq. (|41|) . At A, B chosen according to 



it will be simpler 



/-, x d 2 g 



2K + - - (2K + 2L + l)z 



dg_ 

dz 



+ 



1 



1 



-AM--) -(K + LY 



9 = 0. 



which is of hypergeometric type 

M 1 „ r a M 1 r, r 

a = — -- + K + L, p = -— + - + K + L, 

Again, to get bound states we choose the values 

K = ^>0, L = - € -<0, 

then the quantization rule arises 

a = —N , e N = M + 2N + i> -i. 

It can be noted that €n = e n , when N = n + 1. 

Let us calculate relative coefficient between functions F(z) and G{z). These being 
the form 

F{z) = F z {1+u)/2 (1 - z) {1 -^l 2 F{a, b,c;z), c = | + v , 

M 3 ^ e , M 5 ^ e 
a= T + 4 + 2"2' b= -T + 4 + 2"2 ; 



and 



G(z) = Go z^ 2 (1 - z)^ 2 G(a, (3, r ,z), j = \ + v = c - 1 , 

Mlz/e ^Mlz^e, 
a = T-4 + 2"2 =a - 1 ' ^ = -T + 4 + 2-2 =6 " 1 ' 



must obey the following system 

d 



"" N F+ (+e + M- i/- \)G = , 



/ ; d \/l — Z 



'" X G+(-e + M + v-^)F = 



To find a relative factor, it is convenient to use the second equation 



+ e- 



G 



+{-e + M + u- ^)F = 0. 
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Substituting expressions for F and G, after simple calculation we get to 

dC 1 
2iG — = F (-e + M + u--)F. 

Allowing for the known rule for differentiating hypergeometric functions 

j-0(z) = j- z F(a-l,b-l,c-V,z) 

we obtain 

2*G {a -^\- l) =F Q{ -e + M + ,-\). 



Ultimately, we arrive at the formula 

_ M-1/2 + N 

— * 2 ' 

remembering that eAr = M — 1/2 + 2iV + v. 

5 Radial equations in the case j min 

Let us turn back to the case of the minimal value of j: 
k = +1/2, +1,... 

■j= h - i V*^ h - m h = o , 
-j= h + i h - m h = o ; 



from where for new functions 



H _ h + h q _ h - h 



v/2 ' iy/2 
we derive 

k = +1/2, +1,... 

V$^-H + (-1= + M)G = , 

V¥-?-G - (4= - M)H = . 

And in the same manner for another case we have 
fc = -1/2,-1,... 

/ 4 + i / 4 - M h = , 

^= / 2 - » / 2 - M / 4 = , 
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from whence for new functions 



jj _ h + fi q h — /4 



y/2 ' iy/2 

we obtain 

yf^-^-G + {-A= - M)H = , 

\f$—H - (-A= + M)G = . (56) 
dr 

We can use the above method to eliminate nonphysical singular points. Let us perform 
special transformation on the functions 

G + H = e - p / 2 (g + h) , G-H = e +p/2 (g-h) . (57) 

After simple calculation we arrive at 
instead of (f5lj) 

. d sinh p. . . . , 

— + e— g + ( -e + M - 1/2) = , 
dp cosh p 

d _ e smhp h + (+e + M _ 1/2) (5g) 
ap cosh p 

instead of ([56]) 

. d sinhp, , 

— + e — + -e - M - 1/2) 5 = , 
dp cosn p 

( d_ e ^P ) (+e _ M _ 1/2)/t = . (59) 
ap cosh p 

In the variable z 



r = sinh p 

the system (j58|) takes the form 



, r-Zd e/2 \ , f+e + M-1/2) , , 

Note that the system is symmetric with respect to changes 

f<=>h, e<=>-e. (61) 
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After excluding the function h from (|6U|) we get 

h 2 



(M-l/2) 2 -e 2 ) 



5 = 



(62) 



Ultimately, an equation for g{z) reads 



,<Pg 



2{1 _ z) _| + (1/2 + 



dg /(M-l/2) 2 e 2 + e 1 



5 = 0. 



4 4 1 - z 

In the same manner we get a second order differential equation for h after exclusion of g: 

2 



(63) 
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(+e + M - 1/2) 



(M — 1/2) 2 — e 2 ) 



h = , 



and ultimately 



N d 2 /i . . /(Af-1/2) 2 e 2 -e 1 \, 

Equations ([65]) and ([63]) differ only in the sign at the parameter e. 

6 Solutions of radial equations in the case jmin 

With the use of substitution g = (1 — z) A </?(z), from ([63]) we produce for ip 



z(l - z) + [- - (1 + 2A)z] ^ 



Requiring 



one gets 



^ 



A e 2 + e 



1-z 



A 2 , (M-l/2)- 



2A = e + l,-e 



s // r l / ,n i / 4 -4 2 — (M — 1/2) 2 
z(l - z) if" + [- - (1 + 2,4)z] ^ L_ ^ = o 

w , 1 „ , 4,4 2 - (M-l/2) 2 
</? = r (a, o, c, z) , c = — , a + b = 2A, ab = 



(64) 



(65) 



(66) 



(67) 
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that is 

_ 2A + (M -1/2) 2A — (M — 1/2) 

°~ 2 ' ~~ 2 ' 

Below we will use negative values for A 

A = -e/2, g{z) = {l-z)-^ 2 ^{z)- 

so that 

-e + (M-l/2) , -e-(M-l/2) 

a = , b = . 

2 ' 2 

Any 2-nd order differential equation has two linearly independent solutions; here they 

(pi = Ui(z) = F(a,b,c;z) , 



f2 = U 5 (z) = z l ~ c F{a + 1 - c, b + 1 - c, 2 - c; z) . 
Similar analysis can be performed for eq. (165j) 



^d 2 h . Aln ^dh /(M-l/2) 2 e 2 -e 1 . , , 



With the use of substitution h(z) = (1 — z) L r](z), for rj(z) we produce 



Requiring 



one gets 



z(l-z) rj" + [--(l + 2L)z) r/ 



+ 



L 2 _ £ _ £^-_e\ _±_ _ l2 (M - l/2f 



2 4 J 1-z 



L 



L e 2 -e 



2 4 



2L = +e,-e + 1 



z{l - z) V " + [- - (1 + 2L)z] r, 



, AL 2 - (M -l/2f 



rj = F(a,/3,7,z) , 7 



77 = 
1 



a + /3 = 2L, a/3 



4L 2 - (M-l/2) 5 



that is 



a 



2L — (M — 1/2) 



2L + (M — 1/2) 
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Below we will use negative values for L 

L = (-e + l)/2<0, h(z) = (1 - z)(- €+1 V 2 ri(z) , (76) 

so that 

-e + l + (M-l/2) - e + i-( M -l/2) 

" = g ' " = 2 ' ^ ' 

Functions g(z) and h(z) must obey the above system of first order differential equations. To 
verify that, let us start with the functions 

g = G (l-z) A <p 1 (z), 2A = -e, 
(fl = F(a, b, c,z) c = 1/2 , 
-e + (M-l/2) -e-(M-l/2) 



a 



b= ^ (78) 



h = H (l - z) L m (z) , 2L = -e + l, 
7 ? 2 = z 1 / 2 F(a + l- 7 ,/3 + l-7,2-7,2) , 

a+ l- 7= - t + 2 + f - 1/2) =a + l. 

g +1 . T= -'+»-("-y> =t+1 , 

2 - 7 = c + 1 , (79) 
and relate them with the help of the first equation in the system (|60p 

c2 e/2 \ (-e + M- 1/2) 



After simple calculations we obtain 



G i^-F(a, b, c, z) = ff ( € + M l,2) F{a + 1, b + 1, c + 1, z) , 
az 2 



from whence it follows 



ab (-e + M -1/2) 



c u 2 
that is 

H = i(-e- M + 1/2) G . 
To get polynomial solutions we must require 



note that 



a = -n =>• e n = M + 2n-l/2, 
b = - n -M + 1/2, c = l/2, 

g(z) = (1 - z)-( e " +1 )/ 2 F(a, 6, c, 2) . (80) 



</(*) = (1 - (^-V2)/2 / F( _ nj _ n _ M + 1 j 1 ^) . (81) 
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therefore asz = — r 2 — > — oo the function g(z) tends to zero 

g(z) — > , only if M > - . 

In usual units, that condition for existence of bound states consistent with anti de Sitter 
geometry structure, the inequality M > i, looks as 

p > -— = -A e = 1.213 x 1(T 12 metre 
H 2Mc 2 

so it can be broken only in a very strong anti de Sitter gravitation background, the latter is 
beyond of our treatment. 

Let us write down several energy levels (in usual units) 

„ , 9 1 ch „ , 9 3 ch „ , 9 5 eft 

e = Mc 2 --— , ei=Mc 2 + -— , e 2 = Mc 2 + -— , ... 
2 p 2p 2 p 

(82) 

or 

eo = Mc 2 (l - 1^) , ei = Mc 2 (l + , £2 = Mc 2 (l + ^) , ... 

(83) 

If one mentally increases the curvature radius p, the energy levels will become denser and 
the minimal level tends to the value Mc 2 

e = Mc 2 (l-^) — > Mc 2 . (84) 

7 Conclusions and discussion 

To understand better results, let us discuss the case of minimal j m m in the limit of vanish- 
ing curvature. To this end, let us specify in more detail solutions for minimal values j m m in 
Minkowski space: 



k = +1/2,+!,... 



-1/2,-1,... 



dr 

e fi + i 4- h ~ M h = ; (85) 
dr 



eh + i-4 h-M / 2 = , 
dr 

e h ~ i^r h ~ M h = ° ■ (86) 
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Let detail the case of positive k = +1/2, +1, Let it be 



^ = ft(r) - f -TjT = s{r) (87) 

relevant equations are 

— h+ (e + M) g = , ^-g- (e- M)h = . (88) 
dr dr 

With the substitutions 

h(r) = He^ r , g{r) = Ge 7r (89) 

we get (first let it be (e 2 - M 2 ) > 0) 

7 2 = -(e 2 - M 2 ) == -p 2 , 7 = -ip . 
#7 + (e + M)G = or G 7 - (e - M)# = . (90) 

Thus we have two linearly independent solutions 

h x (r) = H ie +i P r , 91 (r) = G ie + ^ , G 1 = € -^H 1 ; 

ip 

Mr) = H 2 e~ ipr , 92(r) = G 2 e-^ r , G 2 = 6 -^-H 2 . (91) 

ip 

Below, we take Hi = H 2 = 1. We can introduce two linear combinations of these solutions 
the first 

h 1 (r) + h 2 {r) 

= cospr , 

9i{r) + g 2 {r) e-M . 

smpr ; (92) 



the second 



P 



hi(r)-h 2 (r) 

— = smpr , 



2i 

9i(r)-g 2 (r) e-M 



cos pr . (93) 



2i -p 

Now let us specify the case (e 2 — M 2 ) < 0) 

7 2 = -( e 2 -M 2 )== +(7 2 , -y = +q,-q. 
#7 + (e + M)G = or G-/ - (e - M)H = . (94) 

Thus we have two linearly independent solutions 

h^r) = H ie ^ r , 9l (r) = G ie ^ r , G 1 = < —^-H 1 ; 

q 

h 2 (r) = H 2 e-i r , g 2 (r) = G 2 e-i\ G 2 = i^-H 2 . (95) 
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Below, we take H 1 = H 2 = 1. We can introduce two linear combinations of these solutions 
the first 

h 1 {r) + h 2 (r) 

— = cosh qr , 



2 

9i{r) + g 2 {r) e-M 



sinh qr (96) 



the second 



hi{r)-h 2 (r) . 

— = sinh qr , 



2 

9i(r)-g 2 (r) e-M 



cosh qr . (97) 



Evidently, above constructed solutions in de Sitter model provide us with generalizations 
of these of Minkowski space. It may be verified additionally by direct limiting process when 
p — > oo. To this end, let us translate solutions in de Sitter space to usual units 

Ep 

gi (R)=ll + — \ F(a,b,c;—^), c = 1/2 , 

Ep 

f [>2\ 2ch R2 

g 2 (R) = R\1 + F(a + 1 - c, 6 + 1 - c, 2 - c; --^) , 

h 1 (R)=(l + y) F (a,^r,-y), 7 = V2, 

/ jR 2\-H+ 1 /2 R 2 

h 2 {R) = R[l + y\ F(a + 1-7,0 + 1 -7,2 ~T,—^) , 

Parameters of hyper geometric functions are given by 

1 / Ep mcp 1 A 1 / E 1 /? mcp 1 

a= H"^ +( ^"2 ) J' 6= H"^ _( ^"2 



2 V ch h T ) ' 2 V c/i v ft 2 y 

Let us examine the limiting procedure at p — > oo in F(a, b, c; —R 2 /p 2 ) . Because 

1 a(a + 1)6(6 + 1) i? 2 ^4 
2! c(c + 1) ( ~ p 2 ^ ' 
1 a(a + l)(a + 2)6(6 + 1) (6 + 2) i? 2 (p#) 6 

' „2 J 



3! c(c + l)(c + 2) v p 2 ' 6! ' 

and so on, we obtain the following limiting relation 

R 2 

limF(a,6,c; =-) = cospr ==> limg , i(i?) = cospr . 

poo ' poo 
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Similarly, we get 

lim hi (R) = cos pr . (98) 

poo 

In the same manner, we arrive at two limiting relationships 

lim pR Q2{R) = sinpR , lim pR h2(R) = sin pR . (99) 

poo poo 

To rationalize how the finite sums (polynomials of n-order ) may approximate the functions 
cos pR and sin pi? (infinite series), we should take into account the quantization condition 

1 ch 

a = -n =► E = Mc 2 + (2n - -) — 

2 p 

At any fixed E, as p increases the number n also must increase. This means, that the finite 
sums of n terms when p increases will approximate infinite series. 
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